epi draft 



On the competition between the mirror and electromagnetic ion 
cyclotron modes 



R. A. Treumann^'')^'^ and O. D. Constantinescu^ 



003 

O 
O 

(N 

> 
O 

<N 

<D 
O 



Department of Geophysics, Munich University, Theresienstr. ^1, D-80333 Munich, Germany 
Department of Physics and Astronomy, Dartmouth College, Hanover, NH 03755 

Institute of Geophysics and extraterrestrial Physics, Technical University Braunschweig, Braunschweig, Germany 



PACS 96 . 50 . Pw 
PACS 96.50.Vg 
PACS 94 . 05 . Lk 



- Particle acceleration 

■ Energetic particles 

■ Turbulence 



Abstract. - We give a simple argument for the exclusive existence of mirror and elcctromagliotic 
ion cyclotron modes in anisotropic high-/3 plasmas. It is shown that, in addition to a large domain 
of coexistence of both modes, two domains exist in parameter space {A,l3±_) where solely either 
mirror modes or electromagnetic ion cyclotron modes can be excited. In the overlap region the 
modes with the larger growth rate should win. However nonlinear effects may modify such a 
conclusion. 
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Introduction. — The mirror instability is one of the 
* fundamental fluid modes in an anisotropic plasma [2, 9] 
^Siosting an excess of perpendicular energy in the form of 
enhanced perpendicular pressure P±^ > P||. The comple- 
I '"^ nentarv instability is the firehose mode which is excited 
when the above pressure condition is inverted. One con- 
veniently defines a pressure (or temperature) anisotropy 
= (Pj_/f]|) — 1 which quantifies the perpendicular pres- 
p^sure excess. In addition one defines the two components 
, /3_L of the plasma-/?, with (3 = 2fiQNT/ B"^ , as usual. 
Then the simplest necessary linear condition for the mir- 
ror instability [9] to be excited can be written in the ex- 
1-H traordinarily simple form Af3j_ > 1 . More sophisticated 
OO approaches including finite Larmor radius effects [8], fi- 
^^nite electron temperatures [16,17], kinetic corrections and 
^ non-maxwellian distributions [17] and so on, simply mod- 
ify this condition by assigning slightly different mean- 
ings (and in some cases also functional dependencies on 
' other plasma parameters) to the anisotropy A. Surpris- 
ingly, even though the theory of the mirror mode is half a 
century old, its mechanism has still not been satisfactorily 
clarified. One of the reason is that the mirror mode com- 
petes with the electromagnetic ion cyclotron instability 
(EMCI) which grows under exactly the same and appar- 
ently even less restrictive conditions. The necessary con- 
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dition for its excitation is ^4 > 0. One would thus expect 
that the mirror instability looses in growing against the 
EMCI under realistic conditions. 

Observations do however show that the mirror insta- 
bility can exist for a wide range of conditions in many 
places in space plasmas. It has been observed in the mag- 
netosphere, the magnetosheath [24], the solar wind [11], in 
cometary interactions [18] with the solar wind, and in the 
downstream vicinity [5] of Earth's bow shock wave. Most 
intensely it has been investigated using measurements in 
the magnetosheath [4, 7, 14, 21] (for a review of the early 
observations and their interpretation see [19]). Moreover, 
it has been reproduced under specially tailored conditions 
as well in numerical simulations [1, 15, and others]. In the 
observations the mirror mode cannot be detected in its in- 
finitesimal linear state of growth. It is always observed in a 
well developed though evolving time dependent not neces- 
sarily final state when it already has reached large ampli- 
tude and forms deep magnetic holes that are filled with hot 
plasma. Such a state must, however, have emerged from 
the linear state of the instability, and this implies that the 
instability has survived the competition with the EMCI. A 
stationary nonlinear solution has been constructed in [3]. 
Physical arguments for its nonlinear evolution have been 
provided in [12, 17, 23] from different points of view. 

Within the limited range of nonlinear simulations [1,15] 
it has been suggested that mirror modes in the fluid pic- 
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turc form long extended bottles. The stationary analytical 
calculation [3] does not necessarily confirm this conclu- 
sion. The more recent hybrid simulations [1] on the other 
hand show, in agreement with the linear and nonlinear 
theories [13, 17], that mirror modes are short structures. 
Their linear growth rates have been shown to maximize 
in high-/? conditions at transverse wavelengths Aj_ of the 
order of the ion inertial length Ai = c/ujpi, depending how- 
ever on the environmental conditions such that slightly 
longer wavelength may also grow when the short waves 
cannot evolve [17]. Observationally the mirror mode be- 
haves surprisingly. Under very similar conditions in the 
magnetosheath, for instance, it may evolve into short per- 
pendicular and long parallel wavelength magnetic holes 
or also into closely packed magnetic walls. The hybrid 
simulations [1] suggest that small change sin the param- 
eters might produce magnetic holes with no susceptible 
enhancement of the surrounding magnetic field and also 
holes with steep magnetic overshoots at the boundaries. 
Such solutions may apply to the two conditions observed 
in the magnetosheath when the mirror mode evolves. It 
should, however be noted that the magnetosheath plasma 
is bounded from two sides, the bow shock in the upstream 
solar direction and the magnetopause in downstream di- 
rection, and the plasma where the mirror holes are em- 
bedded is flowing. The effect of these boundaries on the 
evolution of the holes is not known and has not been in- 
cluded neither in theory nor in the simulations yet. 

The simpler problem is the competition between the 
EMCI and mirror modes. The competition between the 
two modes poses an unresolved problem which is physi- 
cally minor but requires to be understood. In the magne- 
tosheath both instabilities could grow. In an attempt to 
circumvent the problem [6] the reason for the presence of 
the mirror mode has been attributed to the presence of a 
substantial fraction of He ions, which indeed slightly re- 
duces the growth of the EMCI. On the other hand, the 
mirror mode has also been seen in the absence of He ions. 
However, there is no real need to invoke any complications 
like the presence of heavy ions or strong modification by 
nonlinear theory for finding the parameter ranges of the 
mirror and EMC instabilities. Simple linear theory suffices 
for this purpose. 

Both modes require that the anisotropy A = Pi±/Pi\\ — 
1 > 0. However, the EMCI is a resonant instability 
which depends on the presence of ions with parallel en- 
ergy > miV^/2 = B^/2fioN. This last condition is 
crucial. It has been omitted in all attempts to solve the 
problem where the focus was mostly numerical, consider- 
ing particular parameter combinations in the calculation 
of the growth rates of the two competing instabilities and 
thus restrict to the domain where both are can exist si- 
multaneously. 

For a bi-Maxwellian ion distribution with parallel and 
perpendicular temperatures Tj| , Tl, expressing S^y and the 
number of resonant particles -/Vi.ics through A and i.e. 
eliminating (3n in favour of (3j_ by using the definition of 
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Fig. 1: The existence ranges and thresholds of the EMCI (red) and 
mirror (blue) instabilities in the (A, /3x)-plane. The EMCI grows 
below the red curve, the mirror mode grows above the blue curve. 
At small anisotropics the mirror mode is stable and only the EMCI 
can be excited. In the middle range both instabilities coexist, but 
the one with the larger growth rate will win as suggested in [6]. The 
mirror instability clearly dominates at large anisotropics where the 
EMCI turns out to be stable. 



A, yields that 
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where the right part of the inequality is just the energy 
condition on the EMCI rewritten in terms of A and P±. 
One first observes that, for retaining positive resonant par- 
ticle energies, one requires A < 18 f3± — 1. Solving for the 
conditions on the right restricts this limit even further to 
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The EMCI thus exists only in the anisotropy range < 
A < 2.62 /3_L - 1 and for > 0.38 > /3|| , which is a conse- 
quence of the requirement that sufficiently many resonant 
particles must be present in order to drive the EMC in- 
stability unstable. The condition on P± is unproblematic 
for the mirror instability as the mirror instability refers to 
high-/3 plasma only. 

The condition for the mirror instability, which is an in- 
stability that resonates with the abundant low parallel en- 
ergy (v|| ~ 0) ions [17,22], is simply Ap± > 1, subject 
to the limitations imposed above. The boundaries of the 
regions in {A, P±)-space where these conditions are sat- 
isfied for the mirror and EMC instabilities are shown in 
Figure 1 . There is a large domain between the two curves 



P-2 



Mirror mode 




Fig. 2: The bow shock-magnetosheath system. Shown is the flow 
direction and the evolution of downstream mostly mirror mode tur- 
bulence in the magnetosheath behind the bow shock. The important 
information in this figure aside from the geometry of the structure is 
that the flow is mostly along the magnetosheath thus increasing the 
e-folding time of the mirror instability. The inserts in the foreshock 
show the reflected ions in phase space which when passing ulitmately 
across the shock to downstream both contribute to anisotropy and 
to excitation of EMCI waves just behind the shock, while the fluid 
anisotropy that drives the mirror mode has flrst to evolve along the 
flow. 

where both instabilities coexist and where the one with the 
larger growth rate will win as was suggested in [6]. How- 
ever, at large anisotropies above the EMCI curve solely the 
mirror instability can be excited. With increasing P± the 
anisotropy must become large as well in order still to ex- 
cite the mirror mode. This is clear from the condition that 
the number of high energy resonant particles grows when 
/3_L increases. On the other hand, at small anisotropies the 
mirror instability is stable while the EMCI can grow for 
all /3_L. 

Anisotropies of the order oi A ~ 10 are not unusual 
in the magnetosheath in which case, for moderate P± , the 
external conditions may well lie in the range of suppressed 
EMCI thus favouring the evolution of the mirror instabil- 
ity even in the absence of a substantial fraction of heavy 
ions. Kinetic [17] and nonlinear effects do indeed slightly 
modify the marginal stability boundary and growth rates 
of the mirror instability. In general they lift the lower 
bound of the entire mirror unstable region in Figure 1 up- 
ward displacing parts of this curve into the domain that 
is shared by the mirror and EMC instabilites and thus 
decreasing the unstable mirror range at low anisotropies 
A < 1. The upper mirror unstable range remains less af- 
fected (or even unaffected) by these effects. In that range, 
even for low growth rates 7 <^ LUd of the mirror instabil- 



ity one may expect that the mirror mode can reach fairly 
large amplitudes in a typical convection time across, say, 
the magnetosheath. 

For an estimate we consider magnetosheath conditions. 
The angular ion cyclotron frequency in the magnetosheath 
is of the order of ujd ^ 10 s^^. With a realistic assump- 
tion on 7 ~ O.OlWci, flow velocity of y ~ lOOkm/s, and 
the width of the magnetosheath being of the order of 
~ 2 Re > 10'* km (see Figure 2), a transverse convec- 
tive crossing of the magnetosheath just corresponds to 
one e-folding time of 100 s for a magnetic mirror hole, 
if the hole would form right behind the bow shock. How- 
ever, the anisotropy must first evolve from the bow shock 
along the convective path to reach into the mirror un- 
stable domain. Behind the shock one has (3± ^ 1 — 2. 
Hence, an anisotropy of A '--^ 5 — 10 is needed for the mir- 
ror mode to not compete with the EMCI. The latter we 
know is excited right behind the shock [20]. Anisotropies 
of this magnitude evolve roughly at about half the dis- 
tance between the shock and the magnetopause. Interest- 
ingly enough, it has been found in the observations that at 
about this distance the state of low frequency wave exci- 
tation in the magnetosheath changes about abruptly [10], 
indeed. This apparently reduces the accessible e-folding 
time. However, the flow is not crossing the magnetosheath 
straight from the shock to hit the magnetopause. At the 
contrary, it is deflected and is turned around the magne- 
topause, and observationally the strongest mirror activity 
is found close to the magnetopause at either the flanks or 
at higher latitudes. Thus the flow moves roughly a dis- 
tance of 5-8 Re 6 X 10* km until the mirror modes have 
evolved to the observed amplitudes. This yields a growth 
time of T ^ 600 s or six e-folding times corresponding to a 
factor e^ ^ 400 in the magnetic mirror amplitude. The ob- 
servations show that mirror modes in their evolved state, 
when they are observed in the magnetosheath, have hole 
amplitudes SB— 50-80% (less than) the surrounding ambi- 
ent field value. The latter is of the order of B 30 nT, 
which suggests that the initial infinitesimal mirror wave 
amplitude at about 1 Re behind the bow shock was of the 
order of \SB\ ^ 0.1 — 0.3 nT, which at these low mirror 
frequencies u; ~ is in the measurement noise. This esti- 
mate is consistent with our assumptions on the growth of 
mirrors. 

Nevertheless, the nonlinear state of the mirror mode re- 
mains to be an unresolved problem, as is the basic physics 
of the mirror process. We have elucidated the related prob- 
lems in an earlier paper [23]. In sharpening the argument 
that has been presented there, we note that, in contrast 
to most other plasma instabilities, the mirror mode is an 
exception as the mirrors are no real waves. This contrasts 
for instance its complement, the firehose mode which in 
an anisotropic plasma is a travelling Alfven wave, for all 
its amplitudes, even the largest ones. The mirror mode in- 
stead is a plasma structure. In its infinitesimal initial state 
this does not pose a problem, but once it starts growing it 
is in conflict with thermodynamics as it by itself creates 
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structure, forming a series of magnetic holes in the hot 
plasma with the holes being of the scale of the ion inertial 
length containing grossly unmagnetised ions and magne- 
tised electrons. Structure formation is, however, related to 
phase transitions. Hence the formation of mirror holes in 
hot coUisionlcss plasma resembles an internal phase transi- 
tion of which the plasma is capable. Indeed, the condition 
for instability can be rewritten in the form 

< 2fi„NTj^A EE Bl,, (3) 

where we introduced a critical magnetic field Sera = 
^^2iiqNT±. If the magnetic field in an anisotropic plasma 
drops below this critical value, phase transition sets on, 
and the plasma starts developing structure, breaking off 
into magnetic hole which arc surrounded by regions of 
'normal' plasma state. This argument resembles the Meiss- 
ner effect in superconductivity while it here is applied to 
the ideally conducting (coUisionless) state of a high tem- 
perature plasma. In the same terminology the formation of 
structure on the scale A, <C L, where L is the macroscale 
of the plasma, then means that this phase transition is 
similar to a phase transition of a superconductor of the 
second kind. Of course, this is just a similarity or an anal- 
ogy and not the same physics, as the mirror mode is not 
in a quantum state. It shows, however, that the physics 
of the mirror mode in high temperature plasma contains 
some most interesting problems which are still badly un- 
derstood. 
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